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Abstract. We present analytical and numerical solutions which describe a hot, vis- 
cous, two-temperature accretion flow onto a neutron star or any other compact star 
with a surface. We assume Coulomb coupling between the protons and electrons, and 
free-free cooling from the electrons. Outside a thin boundary layer, where the accre- 
tion flow meets the star, we show that there is an extended settling region which is 
well-described by two self-similar solutions: (1) a two-temperature solution which is 
valid in an inner zone r < 10^'^ (r is in Schwarzchild units), and (2) a one-temperature 
solution at larger radii. In both zones, p oc r~^, Q cx r~^/^, v oc Tp oc r~^; in the 
two-temperature zone, oc r~^/^. The luminosity of the settling zone arises from the 
rotational energy of the star as the star is braked by viscosity; hence the luminosity 
is independent of M. The settling solution is convectively and viscously stable and 
is unlikely to have strong winds or outflows. The flow is thermally unstable, but the 
instability may be stabilized by thermal conduction. The settling solution described 
here is not advection-dominated, and is thus different from the self-similar ADAF found 
around black holes. When the spin of the star is small enough, however, the present 
solution transforms smoothly to a (settling) ADAF. 



INTRODUCTION 



At mass accretion rates less than a few per cent of the Eddington rate, black 
holes (BHs) and neutron stars (NSs) are believed to accrete via a hot, two- 
temperature, radiatively inefficient, quasi-spherical, advection-dominated accretion 
flow, or ADAF [1,2]. While the properties of BH ADAFs are quite well known, hot 
flows onto NSs have not been investigated. Their properties, such as the luminosity, 
spectra, torque applied to a central object, etc., are expected to be different from 
the BH ADAFs because a NS has a surface while a BH has an event horizon [2,3]. 
Here we discuss the structure of a hot accretion flow around a NS [4]. We do not 
attempt a detailed analysis of the boundary layer region near the NS surface. 



SELF-SIMILAR SETTLING SOLUTION 



We consider a steady, rotating, axisymmetric, quasi-spherical, two-temperature 
accretion flow onto a star with a surface, and we use the height-integrated form 
of the viscous hydrodynamic equations. We assume the Shakura-Sunyaev-type 
viscosity parametrized by dimensionless a. We assume viscous heating of protons, 
Bremsstrahlung cooling of electrons and Coulomb energy transfer from the protons 
to the electrons. We neglect thermal conductivity and Comptonization. In the 
inner zone r < 10^'^ (r is in Schwarzchild units, Rs = 2GM/c'^), the flow is two- 
temperature with the density, proton and electron temperatures, angular and radial 
velocities scalings as 

p = Por~\ Tp = Tpor-\ T, = T,or-^'\ Q = Qor-^/\ v = vor\ (1) 

where po? ^po? T^o, fio? are functions of M, a and the star spin s = Q^/Qk{R^), 
and Qk{R) = {GM / R^y/'^ is the Keplerian angular velocity. In the outer zone r > 
10^'^, we have Te = Tp oc and the same other scalings. This self-similar solution 
is valid for the part of the flow below the radius related to the mass accretion 
rate m (in Eddington units, MEdd = 1-4 x lO^^m g/s, and here m = M/Mq): 

m < 2.2 X 10-'al,sl,r;}/\ (2) 

where r^^s = Vs/lO^, ao.i = a/0.1, etc.. The numerical solution of the hydrodynamic 
equations with appropriate inner and outer boundary conditions is represented in 
Figure 1. It is in excellent agreement with the self-similar solition (1). 



PROPERTIES OF THE SELF-SIMILAR SOLUTION 

Spin- Up/Spin- Down of the Neutron Star — The angular momentum flux in the 
flow, J, is negative which implies that the accretion flow removes angular momen- 




FIGURE 1. The radial profiles of density {solid curve), proton and electron temperatures (two 
dotted curves), angular velosity {dashed curve), and radial velosity {long-dashed curve). 



turn from the star and spins it down. This behavior is quite different from that 
seen in thin disks [5,6], where for most choices of the stellar spin parameter s, 
the accretion disk spins up the star with a torque Jthin ~ MQk{R*)RI- In con- 
trast, for the self-similar solution derived here, the torque is negative for nearly 
all values of s. Moreover, J is independent of M. Equivalently, the dimension- 
less torque, j — J/ Jthin, which is ~ 1 under most conditions for a thin disk, here 
takes on the value j ~ — 43m~^Q;^s^ (1 — s'^)^^^ . This torque spins down the NS as 
s — sq/^I + t/r, where the spin-down time is 

T ~ 2 X lO's^jaol {Rm/R,)-^'^ yr or P./P^ ^ 2.7 x \^-^''m^\al^sl^ s-^ (3) 

which is in excellent agrement with observational spin-down rates of some X-ray 
pulsars [7] (here Rm is the magnetospheric radius). Note, the spin-down rate is 
independent of M\ 

Luminosity and Spectrum — The total luminosity has two contributions: from 
the settling flow and from the boundary layer: 

Lss — 6.2 X 10^''mr^^m_2So.i + 8.9 x lO^^mrg ^aQ ^Sg ^ ergs/s, 
Lbl — 1-7 x l{f''^mr^^m^2 ergs/s. (4) 
Note that for sufficiently low m, the luminosity is independent of rh and is dom- 
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inated by the settling flow. Below the radius ~ 45q;o.i Sq.i Comptonization is 
significant (although optical depth is always smaller than unity); therefore the self- 
similar solution (1) is not accurate. The observed spectrum from the settling fiow 
(assuming free-free emission) is calculated to be: 

vL^ ~ 1.7 X 10^^mQ;^.i4i ihv/[2, keV]) ergs/s. (5) 

which is sufficiently accurate up to hv ~ kTg{rc) ~ 400aQ j^^'sq,!^^ keV. 

Stability of the flow — We demonstrate that the settling flow is convectively sta- 
ble and may not have strong winds and outflows (the Bernoulli number is negative) 
if the adiabatic index satisfies 

7 > (3/2) (l - 3^2) / (l - s^A) ~ 1.5. (6) 

The flow is thermally unstable because it is cooling-dominated by free-free emission. 
Stabilization by thermal conduction is studied elsewhere. 
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